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The motivation for the development of a local finite density theory is discussed. One of the 
problems related to an instability in the baryon number fluctuation of the chiral symmetry breaking 
phase of the quark system in the local theory is shown to exist. Such an instability problem is 
removed by taking into account the statistical blocking effects for the quark propagator, which 
depends on a macroscopic statistical blocking parameter e. This new frame work is then applied to 
study color superconducting phase of the light quark system. 



I. INTRODUCTION 

The problems encountered in relativistic high energy physics and astrophysics are local problems about causally 
unconnected parts (in a classical relativistic sense) of the system rather than about the whole system. In both cases, 
the resolution of probing system is much smaller than the size of the systems being observed. The overall feature of the 
system at any given time (in any specific reference frame) is generated by first synchronizing clocks, next measuring 
wanted physical quantity within each part of the system independently at the time that is considered the same within 
the allowed accuracy and then reassemble the whole picture of the system using the bits of local information about 
the system measured (see Fig. |l|). 
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FIG. 1. The process of local measurement in the relativistic space-time. In the first and last steps, the observers remain 
causally (in a classical sense) connected. The observers are separated by space-like distances in the second step when the 
"equal time" measurements are carried out. The correlation between these observers in the second step can only be quantum 
mechanical, which is known to violates classical causality manifested in the violation of Bell's inequality. 



There is no guarantee in quantum world that the refined image of the system obtained under such a procedure 
should be identical to that of the one obtained by doing a low resolution global measurement, in which a uniform (in 
space-time) external field is exerted on the system and the response of the system to the external field is measured, 
like the ones that determine the bulk properties of the system frequently studied in non-relativistic condensed matter 
system. Such a subtle difference is believed to originate from the quantum interference effects which are best illustrated 
by the classical double slit experiments where when one does not care which slit the microscopic particle (e.g., an 
electron) go, there is a interference pattern on the screen behind the double slits, but as long as one determines which 
slit the microscopic particle went through, the interference pattern disappears. The interference pattern can not be 
obtained by summing over the two pictures obtained separately by knowing exactly which slit the microscopic particle 
go. Such a behavior is no longer a theoretical conjecture but a experimental fact which is demonstrated recently by 
experiments from different groups around the world. The same difference should manifest in quantum field theory in 
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the local versus global observables. The importance of an understanding of such a difference has not been appreciated 
due to the lack of a comparing theoretical frame-work. A consistent local theory for quantum field theory at finite 
field strength (leading to finite density of matter and energy) is constructed not only for understand such a 
difference but also for describing the problems raised in high energy physics and astrophysics. 

I shall discuss one of the problems in constructing a local finite density quantum field theory for fermions in this 
talk, those who are interested in the original works can read the papers The particular problem that I am going 

to to discuss here is about the necessity of including fermionic statistical block effects in the chiral symmetry breaking 
phase of light quark system. But before we can discuss the problem, let us first brief review what the local theory 
should be look like at the formal level and what is its most transparent properties. 



II. THE STATISTICAL GAUGE FIELD 

The non-relativistic many body theory for a system with variable particle number is based upon the Grand Canonical 
Ensemble (GCE) The partition function for the system in the GCE is related to its Hamiltonian H and its particle 
number N in the following way 

Z= lim Tre- &-^\ (1) 
n — >oo 

where f2 is the volume of the system, /? = 1/T with T the temperature and fi is the chemical potential. The limit 
£1 — > oo is the thermodynamical limit in which the equilibrium statistical mechanics describes physical many body 
system at equilibrium. A measurement of the ground state fermion number density in GCE which is given by 

Ti =h\l»{l XnZ ) (2) 

corresponds to a global measurement since the chemical potential is independent of both space and time. Global 
measurements are most frequently carried out in condensed matter systems with a size smaller than or comparable to 
our human body (of order 1 meter) which permit us to study the system "as a whole" . The GCE is highly successful 
in describing these situations. Questions posed by high-energy physics, astrophysics and cosmology are somewhat 
different. The size of the systems compared to that of the resolution of the measuring apparatus are normally much 
larger, therefore one actually is capable of study either more details of the probed systems due to the reduction of the 
wave length of the probing system or only small fraction of the probed system due to the impossibility of covering 
the whole probed system by the probing system at the same time. One is in effect doing local measurements in these 
later situations. 

Local measurements in quantum field theory are realized by exerting an external local field to the system at the 
space-time points interested, and, the measured quantity is deduced from the response of the system to the external 
field. For finite density systems, a Lorentz 4- vector local field n a (x), called the primary statistical gauge field, 
is introduced. The ground state expectation value of fermion number density can be expressed as the functional 
derivative 

_. . 5\n.Z 

with Z a functional of [i a . It corresponds to a local measurement in the ground state at the space-time point x. 

For a concrete discussion of the consequences of substituting the global chemical potential in non-relativistic many 
body theory for the local primary statistical gauge field, let us consider an Lagrangian of the following form 

C! = ly(ip + ^0 3 -m)^ + C B + Ant (4) 

with m the fermion mass, Cb the Lagrangian density of boson fields of the system and Ant the interaction between 
the boson and fermion fields. An 8 component "real" fermion field '5 is used for the discussion B with O3 the third 
Pauli matrices acting on the upper and lower 4-component of '5. In most of the cases, Ant = — \&(E(/) — m)^ with 
E a function(al) of the boson fields represented by /. It is well known that the partition functional for the system 
can be formally written in terms of a path integration over the dynamical fields 

Z\n] = J D[f,y,y]exp(i J d 4 xC) = J D[f]exp{iS eS ) (5) 
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which is a functional of p a and lnZ[/i] is its effective action. 

The fermion degrees of freedom can be integrated out first leading to an effective Euclidean action for the boson 
fields and the primary statistical gauge field 

S cS [f, M] = ^SpLn [f f + ^~ S( ? ] + / d±x [C B {f) + up - e] (6) 



2 1 [ip+^0 3 -m] 

with "Sp" the functional trace, e the average energy density of the corresponding free system of fermions (S = 0) and 
p is the ground state expectation value of fermion number density to be discussed in the following at the "chemical 
potential" p = \J~p a ~p a , where ~p a is the ground state value of p a normally provided by the external conditions that 
fixes the density. 

The first term in Eq. ^ is invariant under the U(l) statistical gauge transformation p a (x) — > p' a (x) — p a (x)—d a A(x). 
This is because such a gauge transformation, which transforms the partition functional Z[p] to 

Z[p] Z'[(j] = Z[p'} (7) 

can always be compensated, at least for the first term in Eq. ^ by a change of the fermionic integration variable 

¥(s) -► e iA ^° 3 ^(x) (8) 

*(as) -> e- tA(:E)03 ¥(x) (9) 

so that 

Z\p!\ = Z\ i i\ (10) 

which means that fermionic quantum fluctuation contribution to Z\p] is invariant under the gauge transformation. 
Such a gauge invariance is connected to the conservation of fermion number. Therefore 1 call the 4-vector p a the 
statistical gauge field. Since p a is a local field, its excitation represents certain collective excitations of the system. 
Its general dynamics in different phases of the system at low energy can be derived M . 



III. THE INSTABILITY OF GCE IN THE RELATIVISTIC SPACE-TIME 

The Euclidean effective action for the boson field given by Eq. || can be evaluated in the usual way |Q . Since the 
effective action given by Eq. is a canonical functional of fi a , we can make a Legendre transformation of it, namely 

S e g = S e g - J d 4 xpZ a j a (11) 

with j a the ground state current of fermion number density, to make it a canonical functional of the fermion density 
in order to study the stability of the vacuum state against fluctuations in fermion number. 

The effective potential characterizing the above mentioned energy density can be defined using S e s for space-time 
independent background / fields: 

V e g = -S e g/V 4 (12) 

with V4 the volume of the space-time box that contains the system. The stability of the vacuum state against the 
fluctuations in fermion number density around 7j = can be studied using V e g, which is a canonical function of p. 
Minimization of V e g with respect to p, gives the vacuum value p vac of the system since the vacuum p is a monotonic 
function of p. 

The local finite density theory can be examined by applying it to the chiral symmetry breaking phase of the half 
bosonized Nambu Jona-Lasinio model for the 3+1 dimension and the chiral Gross-Neveu model for 2+1 dimensions 
with its effective potential given by 



, d D p 

V c g = iN q 1 



9 I {2n) D 
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where "D" is the space-time dimension, Go is the coupling constant and p± = (jjq ± p) 2 — p 2 
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The integration over the momentum p^ in Eq. |l3j is divergent, a cutoff is needed to regulate the theory so that finite 
results can be obtained. Whatever the cutoff scheme is, one qualitative feature is true, namely, the stable minima of 
V e ff move from a = to a — ±<r vac ^ as the coupling constant Go is increased above a critical value Gq c . 

Since the logarithmic functions in Eq. ^ has branch cuts on the real p° axis, the choice for the integration contour 
for p° is important. The physical p° integration contour is shown by Cph ys in Fig. ^], which goes under the real axis 
when p° < and above the real axis when p° > 0. 
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FIG. 2. The set of contours belonging to the same physical class. Contour Cph ys is the original physical contour in the 
Minkowski space. Contour Ce is the Euclidean contour. Contour C+ is the quasiparticle contour. ±i7r denote the imaginary 
part of the integrand along the edges of its cuts (thick lines) on the physical po plane. 

It is possible to do the momentum space integration by doing the p° integration first without a cutoff in it, namely, 
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(14) 



In such a case, the integration contours Cph ys , G+ and Ce are all equivalent in the 8-component "real" representation 
for the fermion fields, since the large po behavior of the imaginary part of the logarithmic functions in Eq. ^ is of 
order 0([i 2 /\po\ 2 ) on the physical complex po sheet, which allows the distortion of the integration contour from Cph ys 
to the Euclidean contour Ce and to the quasiparticle contour C+. Let us make a few comment: 1) the 8-component 
"real" representation for the fermion fields make it equivalent between contour Cph ys , G+ and Ce- The usual 4- 
component representation does not has such a property [ fL2| without an arbitrary subtraction. It is in this sense 
that the 8-component theory is not equivalent to the 4-component theory if finite density situation is considered 2) 
the quasiparticle contour G+ is the most convenient contour to demonstrate that the quantum fluctuation term (the 
fermionic determinant) at the right hand side of Eq. |lj is a sum over the quasiparticle energies levels in the negative 
energy sea. The behavior of the quasi-particles is, at least in the approximation taken here, identical to that of 
genuine particles with a global mass despite the fact that the "mass" of the quasi-particles is related to the vacuum 
expectation value of a local field rather than a global constant 3) there is no instability due the fermion number 
fluctuation if the p° integration is done first since only the quasiparticle contributes. This integration contour for 
p° is implied in most of the phenomenological applications in the literature since it is intuitively acceptable. Such a 
"naturalness" is based on an implicit assumption, namely, the quasiparticles is the whole story. 

However, there are a few problems in doing a unbounded p° integration first in a relativistic space-time. The most 
obvious one is related to the fact that a unbounded p° integration, which is done first, corresponds to an arbitrary 
choice of a equal-time hypersurface in space-time with zero thickness. This is not a covariant choice since the spatial 
component of the momentum p^ integration has to be cutted off so that the spatial resolution of the theory is not 
infinitely high. Such a procedure, which is Newtonian in character based on a universal absolute time, is not considered 
fair in relativity since space and time can be related to each other by Lorentz transformation. The second one is 
connected to the observation that the quasiparticle does not saturate the local observables |l],[l3| due to the existence 
of the transient and local quantum fluctuations of the order parameter a. These important kind of transient and local 
quantum fluctuations, which does not decouple in the thermodynamical limit |l3| , are quenched if the momentum p M 
integration is carried out by doing an unbounded p° integration first. 

The relativistic covariance can be maintained by cutting off the theory covariantly in the Euclidean space for the 
momentum p M by restricting the length of the Euclidean p M to be less than or equal to a cutoff A ~ 1 GeV. Such 
a covariance in cutoff is also adopted in lattice simulations of QCD on a symmetric lattice. The covariant cutoff 
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of the theory allow the observers on each space-time point of the measuring hypersurface to have a resolution in 
time (speed of light c = 1) of the same order as the spatial resolution. Therefore we are not talking about absolute 
equal-time, which does not make sense in relativity, but the kind of "equal-time" defined operationally in Fig. |l|. It is 
not expected that the difference between these two kind of cutoff procedure can be removed by renormalization since 
such differences are physical that shall not be renormalized away. 

It can be shown that if the covariant integration procedure is adopted, then d 2 V c tf/dfi 2 \ fl= o < as long as a ^ 
and D > 2, which implies that the a ^ state with ~p ~ is not stable against fermion number fluctuations. Such 
a conclusion is independent of any cutoff scheme (smooth or sharp) and cutoff value for A. The dependence of the 
effective potential on the time component of the statistical gauge field is plotted in Fig. || 
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FIG. 3. The u dependence of V e ff of the conventional approach in the a phase. It has two minima at non-zero u. Here A 
is the Euclidean momentum cutoff that defines the model. 



Such a result is both theoretically and physically unacceptable since it would result in large CP violation in the 
strong interaction vacuum state which is not observed. 



IV. THE STATISTICAL BLOCKING EFFECT 



It appears that the phase in which a vac ^ has a tendency to generate finite density of fermion. Such an unexpected 
behavior may reflect the fact that some thing is missing in our formalism used so far. 

An inspection of the Folk space structure of the vacuum state in the a vac ^ phase reveals the problem. It is know 
that in such a case, the vacuum state is condensed with finite density of fermion-antifermion pairs, namely, 

™QQ = J im -f7^° ( 15 ) 
with the number of the QQ pairs and V3 the spatial volume of the system respectively. This is because 

(?vac ~ (vac \iptp\ vac) (16) 



so 

I vac) = C \ 0) + Ci I QQ) + C 2 \ QQQQ) + ... (17) 

with the number of QQ pairs proportional to V3 if cr vac ^ 0. Like the Fermi sea in the condensed matter system 
of fermions, these fermion-antifermion pairs can produce the statistical blocking effects inside of the system. Such 
statistical blocking effects can not be progressively generated by including higher order loops in a perturbative ex- 
pansion of the effective action. This is because the perturbative terms correspond to a modification of the properties 
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of the system due to finite number of fermions at any specific order of the perturbation expansion. An effect of Eq. 
|j~5| can only be accounted for by non-perturbative means. I proposed a modification of the causal structure of the 
theory to include the statistical blocking effect M in a similar way as the chemical potential do to the theory for the 
finite density situations. Fig. ^.b is a graphical representation of the non-perturbative sea used for our asymptotic 
ensemble in which the maximum filled energy level in the positive energy states is controlled by the statistical block 
parameter e. 
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(a) The Perturbative Dirac Sea 



(b) The Non-perturbative Dirac Sea 



FIG. 4. The perturbative and non-perturbative Dirac sea for the bare particles. The hollow circles represent the empty 
states. 

It can be shown that if Fig. ^.b is chosen, the energy integration contour Cp} lys has to be distorted to Cph ys drawn 
in Fig. | 
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FIG. 5. The p° integration contour for the effective action for transition amplitudes between states in which both the 
fermion and antifermion states with absolute value of their energy below e filled. The thick lines extending to positive and 
negative infinity represent the branch cuts of the logarithmic function. 

The effective potential 14// depends, in addition to other conventional ones, the statistical blocking parameter e. 
Eq. O should be replaced by 



V e s(a, fj,, e) = iN g 



d D p 



In I 1 - 
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(18) 



which should be evaluated in the Euclidean space in a way that preserves the analytic properties of Fig. |5| and cutted 
off covariantly. 
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The stability of the resulting effective potential for D = 4 is represented in terms of equal V e ff contour in Fig. g 
with a nonvanishing a/A = 1 and A ~ 1 GeV is the covariant cutoff. The effective potential is not stable if both /i 
and e are zero. The minima at finite fi are now local minima rather than global. The global minima are located at 
finite e and vanish fi at which point d 2 V e s/dfj, 2 > and d 2 V c s/de 2 > which means that they are stable. 



0.15 
0.1 
0.05 

o [L/A 

-0.05 
-0.1 
-0.15 

-0.15 -0.1 -0.05 0.05 0.1 0.15 

8/A 



FIG. 6. The contour plot of the equal V e g in the phase in which D — 4, a/A = 1.0. The point /i = and e = is not stable 
both in the ±e direction and in the ±/i direction. The absolute minima of the effective potential is located at about /i = and 
e/A = ±0.1. There is also a pair of local minima located at fj,/A = ±0.1 and s = 0. 

Therefore, the statistical blocking parameter cure the CP problem mentioned above. Some of the implications of 
having a non- vanishing e in the a vac 7^ phase of the vacuum is discussed in Refs. |]||. 

V. THE PROPERTIES OF COLOR SUPERCONDUCTIVITY 

The ground states of the strong interaction could have different phases from the chiral symmetry breaking phase 
or the a-phase 0-^]. They are characterized by a condensation of diquarks causing color superconductivity. Such a 
possibility is interesting because it may be realized in the early universe, in astronomical objects and events, in heavy 
ion collisions, inside nuclcons [f7|,p|jic|l and nuclei, etc.. For the possible scalar diquark condensation in the vacuum, a 
half bosonized model Lagrangian is introduced 0,0 , which reads 

Ci = ^[ip-a-in- f 7 5 3 - l 5 A cX c O {+) - 7 5 ^ c X c O(_)] * - -±-(a 2 ± tt 2 ) ± ^;X C X C , (19) 

where a, if, x c and \ c are auxiliary fields with (x c )^ = —X c an( i ^Oj G3/ are coupling constants of the model. A c and 
A c (c = 1,2,3) act on the color space of the quark; they are A c ClC2 — —e CClC2 A C:Cl c 2 = e CClC2 with e abc (a,b = 1,2,3) 
the total antisymmetric Levi-Civita tensor. Here 0(±) are raising and lowering operators respectively in the upper 
and lowering 4 components of 

This model has two non-trivial phases. The vacuum expectation of a is non-vanishing with vanishing \ 2 = ~x c X c 
in the a-phase. The vacuum state in the a-phase is condensed with quark-antiquark pairs. The vacuum expectation 
of a is zero with finite \ 2 that spontaneously breaks the U(l) statistical gauge symmetry in the second phase, which 
is called the w-phase. There is a condensation of correlated scalar diquarks and antidiquarks in the color 3 and 3 
states in the w-phase of the vacuum state. 

Since diquarks and antidiquarks are condensed in the w-phase, it is expected that an exchange of the role of e and 
(i occurs. It is found to be indeed true: there are also two sets of minima for the effective potential, the first set is 
the one with finite [i = ±[i vac and e = and the second set corresponds to /i = and finite s = ie vac in the w-phase 
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of the model. But here the absolute minima of the system in the w-phase correspond to the first set of solutions, in 
which the CP invariance and baryon number conservation are spontaneously violated due to the presence of a finite 
vacuum [i% ac . The stability of the resulting effective potential is represented in Figs. [?] and || 




e/A 



FIG. 7. The contour plot of the equal V^g- in the phase in which x-uac/A = 1.0. The point fi = and e — is also not stable 
both in the ±s direction and in the ±fx direction. The local minima of the effective potential is located at about fi — and 
e/A — ±0.15. There is also a pair of absolute minima located at fi/A — ±0.3 and e — 0. 




FIG. 8. The 3D plot of V eS in the phase in which Xvac/A = 1.0. 

This conclusion is also applicable to the /3-phase of p-wave diquark condensation models ^|,|| with vector fermion 
pair and antifcrmion pair condensation, in which the chiral symmetry SU(2)l x SU(2)r is also spontaneously broken 
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down. It is shown in Fig. || and |h] 




FIG. 9. The contour plot of the equal V e Q in the phase in which \J 4>% ac / 'A = 1.0. The point /i = and e = is also not 
stable both in the ±e direction and in the ±/Lt direction. The local minima of the effective potential is located at about /i = 
and e/A = ±0.15. There is also a pair of absolute minima located at fi/A = ±0.3 and e = 0. 




FIG. 10. The 3D plot of V eS in the phase in which a/^Lc/A = 1.0. 

Therefore, the statistical blocking parameter e and the primary statistical gauge field /i Q provide a more complete 
set of macroscopic parameters to describe the condensation of fermionic particles in a relativistic system. Their 
physical meaning is also clear. The statistical gauge field characterizes the condensation of particle (pairs), which has 
net fcrmion number. The color superconducting phase is caused by condensation of fermion pairs, which implies that 
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fermion pairing or antifermion pairing are favored against the fermion-antifermion pairing. The matter antimatter 
tends to separate in the color superconducting phase resulting in genesis of matter. The statistical blocking parameter 
characterizes the concentration of fermion-antifermion pairs with no net fermion number. The chiral symmetry 
breaking phase, on the other hand, favors the fermion-antifermion pairing. It is not surprising that e ^ in the chiral 
symmetry breaking phase. 



VI. SUMMARY 



In summary, it is important to take into account, among others ones discussed in Refs. |l|j2|], of the statistical 
blocking effect in a localized finite density theory. 
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